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1 Abstract 



We derive certain constraints on the reflection matrix for reflection from a 
plane, nonmagnetic, optically anisotropic surface using a reciprocity theorem 
stated long ago by van de Hulst in the context of scattering of polarized light. 
The constraints are valid for absorbing and chiral media and can be used as 
tools to check the consistency of derived expressions for such matrices in 
terms of the intrinsic parameters of the reflecting medium as illustrated by 
several examples. 

2 Introduction 

Except in the simplest situations where linearly polarized light is incident 
on an optically isotropic surface or a specially oriented anisotropic surface, 
the polarization state of a light wave in general changes when it is reflected 
off a surface. Changes in polarization states imply changes of phase, as 
we know from the work of Pancharatnam. A correct analysis of optical 
devices involving polarization, with or without interference, requires a precise 
and reliable method for handling polarization transformations resulting from 
reflections. These are described by means of a 2 x 2 complex reflection matrix. 
The specification of this matrix requires the choice of a set of basis states 
for the incident as well as the reflected waves alongwith their phases. In this 
communication we first describe the most natural convention for the basis 
states and a consistent one for the description of reciprocity. We then derive 
certain constraints that reflection amplitudes must satisfy on account of the 
principle of reciprocity and illustrate their use with several examples. The 
principle of reciprocity in polarization optics has been widely discussed in 
literature and is different from time-reversal invariance in that (i) it applies 
to systems with absorption and (ii) it deals with one incident wave and one 
scattered wave at a time. For a recent review we refer the reader to Potton [1] . 
For the purpose of this paper the most appropriate formulation of reciprocity 
is the one given by van de Hulst [2] wherein he choses to rotate the scatterer 
instead of reversing of the wave. 
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3 The phase convention 



It was pointed out recently [3l H] that when a plane wave of light changes 
its direction of propagation from a wave vector to a wave vector k' due to 
refraction, reflection or scattering, the definition of the matrix that relates 
the incident polarization state to the final polarization state, i.e. the Jones 
matrix, requires the choice of a set of basis states and their phases for each of 
the two directions of propagation. The choice that is most often used both 
in polarization optics as well as in scattering theory is the following, also 
illustrated in Figs. 1 and 2. 

For the wave vector k, a set of orthogonal, linearly polarized states along 
X and y, called the p and the s states, are chosen as the base states 1 and 
2 respectively where x is in the plane and y is perpendicular to the plane 
of reflection or scattering. The phases of the basis states are chosen such 
that in the basis state 1, = Eexp{iujt) , Ey = and in the basis state 2 
Ere = 0, Ey = Eexp{iujt), where E^ and Ey are the x and y components of 
the electric field in the wave. With this convention, the vector (l/\/2)col.[l,l] 
represents a linearly polarized state along a direction lying in the (x, y) plane, 
making an angle 45° with x and the vectors (l/A/2)col.[l,±i] represent the 
right and left circularly polarized states. For the wave vector k', the con- 
vention most often used in scattering theory as well as in polarization optics 
is the following: Rotate the basis states about an axis perpendicular to the 
plane of reflection, i.e. along y, through an angle such that k goes to k'. Let 
X and y go to x' and y' under this rotation. The polarization basis states for 
k' are chosen to be linearly polarized states along x' and y' with their relative 
phases chosen such that in the basis state 1, E'^ = Eexp{iut), Ey = and 
in the basis state 2 E'^ = 0, E'^ = Eexp^iut). Note that since the rotation 
is about y, y' = y. This choice leads to a precise phase convention for back- 
scattering or for reflection at normal incidence at a surface, where k'=—k. 
For this case one gets x' = —x and y' = y. We shall call the above convention 
the "travelling-frame convention". In literature the above described choice 
of phase convention is referred to as "choice of the coordinate system" . 



4 The reciprocity constraints 

Scattering of a polarized plane wave with wave-vector k to a wave with 
wave-vector k' from a scatterer is described by a 2 x 2 complex scattering 
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matrix A whose matrix element Aij represents the complex amplitude for 
an incident wave with unit amplitude in polarization state i to scatter into 
the polarization state j. If A{k, k') represents the matrix for scattering from 
k to k' and A{—k', —k) the matrix for the reverse scattering, the correct 
statement of the principle of reciprocity with the above phase convention has 
been given by Sekera [S] as 

A{-k', -k) = A{k, k') (1) 
where the matrix A is the "n-transpose" of A, defined [3] as, 

A, = (-1)^+^,. (2) 

For a 2 X 2 matrix, A is the transpose of A with a change of sign of the 
off-diagonal elements. 

In order to derive the constraints on the reflection matrix resulting from 
reciprocity, a somewhat different formulation of the principle of reciprocity, 
first made by van de Hulst [2] in the context of scattering problems, is more 
useful [6]. van de Hulst choses to rotate the scatterer rather than reverse 

— # — * 

the direction of propagation. Let the polarization basis states for k and k' 
be defined using the travelling-frame convention, let k and k' be along z and 
z' respectively and let the scattering be in the {x,z) plane as shown in Fig. 
1 where a is the scattering angle, van de Hulst's theorem states: If the 
scatterer is rotated through 180° about an axis defined by the line bisecting 
the angle between the vectors k' and —k , called the bisectrix, the matrix 
A{k, k') goes to the matrix A{k, k') where A is the n-transpose of A defined 
by Eq.©. 

The theorem, translated to the problem of reflection from a plane surface 
in optics, can be phrased as follows: If the reflecting medium, assumed to 
be reciprocal, is rotated about the normal h to the surface SS through 180° 
(Fig. 2), the reflection matrix Z goes to Z' where 

Z' = Z, (3) 

where Z is the n-transpose of Z defined by Eq.([2]). 

The theorem is true in the presence of absorption and dichroism and has 
the straightforward consequence that if the reflecting medium, assumed to 
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be reciprocal, is invariant under a rotation through tt about fi, the reflection 
matrix Z, for any angle of incidence, must be antisymmetric. For such cases 
therefore, 

Zij = -Zji. (4) 

Such cases include: 

A. An optically isotropic medium, i.e. a medium with no birefringence or 
dichroism, linear or circular. 

B. A medium with only optical activity and circular dichroism but no 
linear birefringence and no linear dichroism, 

C. An absorbing uniaxial medium with or without optical activity, with 
the optic axes for birefringence and dichroism coinciding and being perpen- 
dicular to the surface, 

D. An absorbing uniaxial medium with or without optical activity, with 
the optic axes for birefringence and dichroism coinciding , and lying in the 
plane of the surface, 

E. A nonabsorbing biaxial medium with one of the principal axes perpen- 
dicular to the surface. 

In addition to the above cases, when light is incident normally, any re- 
flecting medium is invariant under a tt rotation about the direction of the 
incident beam. At normal incidence therefore, the reflection matrix for any 
reciprocal medium must be antisymmetric. 

In cases A, B and C, when light is incident on the surface normally, 
there is an additional constraint when the reflecting surface is invariant under 
rotations about the direction of incidence, i.e about the normal to the surface. 
In the travelling frame convention this additional constraint can be expressed 
as, 

R{^)ZR{4)) = Z, (5) 

where 

^(0)=f''t 1 (6) 

y sm0 COS0 j 

is a rotation matrix that rotates an incident Jones vector by an arbitrary 
angle (j). 
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Equation can be proved as follows: Let the reflection matrix of the 
unrotated sample be Z and let | ipf > be the flnal state resulting from 
reflection of an incident state \ ipi > so that 



(7) 



The reflection matrix Z of the rotated sample is obtained from the 
condition that when the state | tpi >, rotated through an angle 0, is incident 
on the rotated sample, the reflected state, in the travelling frame, must be 
the state \ ipj >, rotated by an angle —(f) i.e. 



Since is arbitrary, requirement of invariance under rotation about the beam 
axis therefore gives Eq. ([5]). 

Let the matrix that satisfles the reciprocity constraint (jlj), as well as the 
isotropy constraint be called Zq. It can easily be shown that Zq must be 
of the form 



where r is a complex number. In other words whenever a reciprocal reflecting 
medium is invariant under an arbitrary rotation about the normal to the 
surface its reflection matrix for normal incidence is given by Zq. The result 
is known to be true for optically isotropic surfaces. The fact that it is true 
in the presence of optical activity and circular dichroism and that it follows 
from simple symmetry considerations came as a news to the author. We also 
note that since the matrix Zq is diagonal it is robust against a phase change 
between the two basis states. 




(8) 



Eqns.© and give 




(9) 




(10) 
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5 Dependence on the polarization basis 



We next discuss the dependence of the theorem given by Eq. ([3]) on the basis 
states used to express the Jones vectors for the propagation directions k and 
k'. Although while stating the theorem we assumed a basis of "in-phase" 
linearly polarized states along x and y, this is by no means the only possible 
choice for the theorem ([3]) to be true. One can also choose as basis states, a 
pair of orthogonal elliptically polarized states with the principal axes of the 
polarization ellipses being along x and y; the states being phased such that 
at t=0, in the wave with wavevector k, the basis state 1 has Ey = and the 
basis state 2 has = 0. Similarly, in the wave with wavevector k', at t=0, 
Ey' = in basis state 1 and = in basis state 2. The theorem as stated 
above remains valid in this set of basis states. This can be proved easily. 

The basis described above is obtained from the original linearly polarized 
basis by means of a unitary transformation U that takes the state | x > 
i.e. the state with coordinates (90°, 0°) on the Poincare sphere, along a 
geodesic arc, to a state | P > with coordinates {90° + 1], 0°), where 90° > r] > 
—90°. Such a transformation is achieved by means of a linearly birefringent 
waveplate with retardation rj, with its fast axis at 45° to the x axis. The 
matrix U is therefore given by 

U = = R{^h)Lo{r])R{-Ah). (11) 

where Lp{ri) is the Jones matrix for a linear retarder with retardation rj and 
fast axis making an angle f3 with x. The fields in the basis states 1 and 2 are 
given by, (1) E^ = Ecos^f] /2)exp{iujt) , Ey = —iEsm{r] / 2)exp{iut) and (2) 
Ere = iEsin{ri/2)exp{iujt), Ey = Ecos{ri/2)exp{iut). The cases t] = and 
7] = tt/2 give the fields in the linearly polarized and the circularly polarized 
basis respectively. 

The matrices Z and Z' when transformed to the new basis are given by 
Q and Q' where 

Q = UZU^ and Q' = UZ'U^ (12) 
It can easily be shown by required matrix multiplication that 

Q' = Q (13) 
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Equation (fT^ states the reciprocity principle in the basis of the chosen elhpti- 
cally polarized states. It needs to be mentioned however that in an elliptically 
polarized basis, the form of the matrix under conditions of normal incidence 
and spatial isotropy ([5]) is not given by Eq. ffTOl) . The latter requires R{4>) to 
be of the form ([6]) which is true only in the linear basis. 

6 Applications 

The constraints derived above can be used as tools to check derived expres- 
sions for the matrices of reflection from optically anisotropic surfaces in terms 
of the intrinsic parameters of the sample. Under appropriate conditions the 
derived expressions must satisfy these constraints. We cite below some exam- 
ples from literature where derived expressions for reflection matrix elements 
indeed do so. 

Sosnowski [8] has derived the reflection matrix elements for reflection 
from the surface of a uniaxially anisotropic medium placed in an isotropic 
ambient medium for the case when the optic axis is parallel to the interface 
and is oriented at an angle a from the plane of incidence. These have been 
reproduced on p. 355 of [7J. First consider the case of normal incidence i.e. 
(/)o = 0. We derived the expressions for the off-diagonal elements Vps and Vgp 
for this case using the formulae in Eqs. (4.244) - (4.246) of [7j. It was found 
that they satisfy r^^ = —rsp as required by Eq.(jl]) of this paper. Next consider 
the case of oblique incidence i.e. 0o 7^ 0. For this case we programmed the 
above chain of formulae on an Excel worksheet and computed Vps and Vsp 
for several hundred randomly chosen sets of the parameters Nq, Nio, Nie, 
a and 0o in their allowed ranges where Nq is the refractive index of the 
isotropic ambient and Nio, Ni^ are the two refractive indices of the anisotropic 
medium. In every case we obtained Vps = —rgp. To cite two specific examples, 
for no = 1.2, riio = 1.7, riie = 1.3, 0o = 30° and a = 60°, we obtained 
Tps = .0637 and r^p = —.0637 and for uq = 1, riio = 1.2, riie = 1.5, (po = 75° 
and a = 10°, we obtained Vps = —.0274 and r^p = .0274. 

Engelsen |9j has derived the expressions for the matrix of reflection from 
a uniaxially anisotropic film on an isotropic substrate in an isotropic ambient 
medium for the case where the optic axis of the uniaxial medium is perpen- 
dicular to its boundaries with the substrate and the ambient. These have 
been reproduced on pages 356-357 of [7]. The matrix is diagonal in this case. 
We derived the expressions for the diagonal elements r^^ and Vpp for normal 
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incidence using the formulae in Eqs. (4.249) - (4.257) of [7]. It was found 
that they satisfy r^p = —r^s as required by Eq. lfTII]) of this paper. 

We next consider some examples from the literature on reflection from 
a reciprocal, isotropic chiral medium where the constraints derived in this 
paper yield useful insights. 

Lekner [TU] has derived expressions for the reflection matrix for reflection 
from the boundary of an achiral and an isotropic chiral medium. For normal 
incidence these expressions yield Vpp = r^s and Vps = Vgp = 0. For oblique 
incidence the expressions satisfy Vps = Tgp. These results differ from the 
results of this paper by a sign. The reason lies in the phase convention for 
the reflected p-wave used in [TU] (see section 1.3) which differs by vr from 
the one used in this paper (Fig. 2), resulting in a change in the sign of the 
reflected p-wave amplitude. This changes the signs of r^p and Vpp. When this 
change of sign is accounted for, the results of ^lOj agree with those of this 
paper. 

Silverman [11] derived the reflection matrix for reflection at the surface of 
an isotropic, nonmagnetic chiral medium for two sets of constitutive relations 
that are (I) invariant and (II) noninvariant under a duality transformation 
of the electromagnetic fields. The symmetric constitutive relations (I) lead 
to null differential reflection at normal incidence for incident right and left- 
circular ly polarized light. The asymmetric constitutive relations (II) on the 
other hand lead to non-zero differential reflection for right and left-circularly 
polarized light. The author indicates a preference for (I) based on some 
difficulties with the results obtained from (II). Using the constraint stated 
above, i.e. the reflection matrix for this case must be given by Eq. ffTOj) . 
any theory that yields non-zero differential reflection at normal incidence 
for incident right and left-circularly polarized light can be ruled out. If 
we assume that the derivations in [TT], which do satisfy our constraints, 
are correct, it could be concluded on grounds of symmetries alone that the 
asymmetric constitutive relations are incorrect. 

Georgieva [12] reported a solution for the amplitudes for reflection from 
the surface of a reciprocal optically active medium using a corrected Berre- 
man's matrix, arguing that Berreman's matrix is incorrect since it yields 
unequal off-diagonal elements for the reflection matrix. The considerations 
of this paper support this assertion. Interestingly however, while the off- 
diagonal elements in [12] are equal and opposite in sign as required by the 
above constraints, the diagonal elements do not satisfy these constraints 
when light is incident normally. Eqns.(27) and (30) in [12] yield, for nor- 
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mal incidence, Vgs = Tpp. The constraint given by Eq. fllOl) however requires 
= —Tpp. The negative sign is non-trival as it represents the difference 
between a plane glass plate and a halfwave plate. In this work since Vss = Tpp 
and Tps = —Tsp, the disagreement cannot be explained by a phase convention. 
We believe therefore that there is a problem with the derivation in [12] . 

7 Discussion 

Our reason for dwelling at length on the conventions regarding basis states is 
that the reflection matrix as well as the statement of the reciprocity principle 
depend on these conventions. While the travelling-wave convention is a fairly 
standard one, used by most textbooks on optics [131 E] to relate the polar- 
ization states for k to those for k' for defining the reflection matrix, there 
are occasional exceptions. For example Lekner [10] and Bassiri et al. [H] 
use a different convention which we shall call the "fixed-frame convention". 
Consequently they obtain expressions for the Fresnel reflection amplitudes 
for reflection off a chiral surface that differ from those in [T3| [7] even in the 
limit when the chiral parameter goes to zero. As mentioned before the am- 
plitudes obtained with the two conventions are related by a change of sign 
of the amplitude of the reflected p-wave, hence of r^p and Vpp. 

In the analysis of propagation problems involving a series of oblique re- 
flections terminating in a reflection at normal incidence so that the beam 
retraces its path, as for example in a Michelson interferometer, the problem 
of phase convention occurs twice, once while defining the reflection matrix 
and again while relating the forward and backward propagating waves. Since 
the first choice implies a choice for normal incidence, the natural thing to 
do is to make the second choice to be consistent with the first one. Unfor- 
tunately this has not always been the practice in literature. For example in 
Vansteenkiste et al.[T5] the travelling frame convention is used for the re- 
flections at oblique incidence and a fixed-frame convention for relating the 
forward and backward propagating waves. As a consequence the matrix for 
reflection at normal incidence is defined differently from the ones for oblique 
incidence. We find this somewhat unsatisfactory and that it is avoidable 
if one consistently uses the travelling-frame convention. As demonstrated in 
[T5] it is indeed possible to derive correct results if one carefully keeps track of 
the phase conventions. However both in regard to pedagogy and applications 
it would be desirable and simpler if a consistent convention were used and all 
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reflections described similarly. If the travelling frame convention is used con- 
sistently the matrices for reverse propagation are of course n-transpose of the 
corresponding matrices for forward propagation instead of being the trans- 
pose |4] . Though less familiar, the n-transpose is, however, an equally simple 
and elegant mathematical construct that satisfies the property {AB) = BA. 

The use of the fixed-frame convention for reflection amplitudes has some- 
times been justified by arguing that for normal reflections off an optically 
isotropic surface it yields a unit reflection matrix which avoids the asymme- 
try between the s and p wave reflection amplitudes. We point out that this 
is achieved at the expense of counter-intuitive behaviour of the amplitudes 
elsewhere. For example for reflection off ideal metallic mirrors at grazing 
incidence the fixed-frame convention gives a unit matrix suggesting no po- 
larization change. We know however that under these conditions a right 
circularly polarized wave is reflected as a left circularly polarized wave and 
vice versa. Another problem with the use of the fixed-frame convention is 
that there is an asymmetry of conventions between the transmitted and the 
reflected waves. In a scattering problem there is no natural place for such 
an asymmetry. The neat correspondence between the theory of scattering of 
polarized waves and that of reflection and refraction is thus needlessly given 
up. 

To sum up, in the examples discussed above we found cases ([8], [9] and 
[llj) where the derived expressions satisfy the constraints derived in this pa- 
per. We found cases ([TU] and [T^) where they do so after accounting for 
a difference in phase convention. Finally we found ([12]) where the 

derived expressions donot satisfy the constraints and we conclude that the 
derivation is incorrect. We wish to emphasize however that the satisfaction 
of the constraints is a necessary but not a sufficient condition for the correct- 
ness of the derived reflection amplitudes. The constraints therefore provide 
only a partial test for the derived amplitudes. Finally we note that all the 
considerations in this paper relate to the linear regime of optics and do not 
include nonlinear phenomena. 
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BISECTRIX 



Figure 1: The geometry of scattering of a plane wave with wave vector k 
along 5 to a wave with wave vector k' along z' where z and z' lie in the X-Z 
plane. The coordinate system {x',y',z') defining the polarization basis states 
in the scattered wave is obtained from the {x,y,z) system in the incident wave 
by a rotation about y through an angle a which is the scattering angle. 
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Figure 2: The geometry of reflection of a plane wave propagating along z from 
a plane surface SS whose normal along h lies in the X-Z plane. The angle of 
incidence is d and the relation between the coordinate systems {x,y,z) and 
{x',y',z') is the same as in the scattering problem illustrated in Fig. 1. 
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